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Abstract 

It is shown in this paper that non-conforming finite elements on the 
triangle nsing P^-nonconforming polynomials and -conforming poly¬ 
nomials can be easily built and used.They appear as an ’enriched’ ver¬ 
sion of the standard piecewise quadratic six-node element.This work is 
divided into two parts.In the first we present the basic- property of the 
element,namely how it can be built and basic error estimates.We have ob¬ 
served that this new element behaves like non-conforming element.In 
the second part we have applied our element to the elliptic problem and 
the theoretical estimate has been guaranteed by numerical result. 


1 INTRODUCTION 

The finite element method has achieved great success in many fields of science 
and technology since it was first suggested in elasticity in the fifth decade of 
20*^ century.Today it has become a powerful tool for solving partial differential 
equations mm .The key issue of the finite element method is using a discrete so¬ 
lution on the finite element space,usually consisting of piecewise polynomials,to 
approximate the exact solution on the given space according to a certain kind of 
variational principle. When the finite element space is a subspace of the solution 
space, the method is called conforming.lt is known that in this case,the finite 
element solution converges to true solution provided the finite element space 
approximates the given space in some sense. 

In general, for a 2m order elliptic boundary-value problem,the conforming 
finite element space is a subspace.lt means that the shape function in 
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this conforming finite element space is continuous together with its m — 1 order 
derivatives.That is for a second-order problem, the shape function is continuous 
and for a fourth order problem, the shape functions and its derivatives are 
continuous.lt is a rather strong restriction put on the shape functions in the 
latter case. 

It is prove that to build up a conforming finite element space with C^- 
continuity for a two dimensional fourth-order problem, like the plate bending 
problem in elasticity, at least a quintic polynomial with 18 parameters is re¬ 
quired for a triangular element and a bi-cubic polynomial with 16 parameters 
for a rectangular element.lt causes some computational difficulties because the 
dimension of the related finite element space is fairly large and its structure is 
rather complicated. 

Therefore, it is desirable to relax directly the C'^~^ continuity of the finite 
element space.lt comes to the so called ’nonconforming’ finite element method 
which had and still has a great impact on the development of finite element 
methods.However, it was found that some nonconforming element elements con¬ 
verge and some do not. The convergence behaviour sometimes depends on the 
mesh configuration. 

Up to now there has been proposed a vast number of engineer devices based 
on different mechanical interpolations,like unconventional elements,energy-orthogonal 
elements with free formulation,quasi-conforming elements,generalised conform¬ 
ing elements and many others.The approximate spaces related to all these ele¬ 
ments mentioned above are not included in the given solution space.Hence they 
are simply called ’nonstandard’ finite elements.A unified mathematical treat¬ 
ment for analysis of these nonstandard finite elements has been proposed by 
various authors. 

A brief description and analysis of some interesting and important conform¬ 
ing and nonconforming finite elements has been given below. 

nonconforming element mm- This is a triangular element which is not 
. The shape function is a linear polynomial with three nodal parameters at 
mid points of three edges of the triangle.This element converges for second order 
elliptic problem with optimal rate. 

Wilson-element^ 01 [S]. This is a nonconforming rectangular element.The 
shape function is a quadratic polynomial with six parameters, four at vertices 
of the rectangle and two internal degrees of freedom, like the second order 
derivatives.This element converges for rectangular mesh, but does not converge 
for arbitrary quadrilaterals.lt is interesting to mention that the behaviour of the 
Wilson element is better than the corresponding bilinear Qi conforming element 
as many engineering example have indicated. 

The rotated Qi element liiTiiH]. This is a newly established non-conforming 
rectangular element.The shape function consists of four terms as[l,x,y,a;^ — 

There are two versions of choosing nodal parameters.The first one uses four 
function values at the mid-point of each edge of the rectangle.The second ver¬ 
sion uses four mean values of the shape function along edges. Both versions are 
convergent for rectangular meshes.However,the first version is not convergent 
for arbitary quadrilaterals unless certain mess conditions are satisfied. 

Morley element [5] This is an old and simplest plate element.The shape func- 
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tion is a quadratic polynomial with six nodal parameters.They are three function 
values at vertices and three normal derivatives at mid-points of three edges.This 
element does not even belong to (7° class,nevertheless it is convergent for the 
fourth-order problem[Tn] . Surprisingly, it is recently proved that the Morley 
element is divergent for the second-order elliptic problem la¬ 
in contrast, it is well known that there exists for long time the conforming 
element for the second order problem.Its shape function is again a quadratic 
polynomial with six parameters as three function values at vertices and three 
function values at mid-points of edges.This quadratic element is convergent 
for the second order problem, but divergent for the fourth-order problem. 

Zienkiewicz incomplete cubic triangular element [12j . The shape function 
consists of incomplete cubic polynomials with specially chosen nine terms.The 
nine nodal parameters are three function values together with six hrst partial 
derivatives with respect to x and y at three vertices.This is a C° element but not 
.It is proved that this element is convergent only for very special meshes, 
namely, all edges of triangles are parallel to the three given directions.lt is 
a very interesting phenomenon that the Zienkiewicz element using the cross¬ 
diagonal mesh actually tend to a limit,but it is not the true solution of the 
given problem,rather of another problem |13j . 

Another new nonconforming piecewise quadratic finite element on triangles 
has been discussed in [14]. This element satisfy patch test of Irons and Razza- 
qque m .This implies that on element interfaces,one should ensure the continu¬ 
ity of the approximation at the Gauss-Legendre quadrature points needed for 
the exact integration of third-degree polynomials. Optimal-error estimates and 
regularity properties for Dirichlet’s problem has been studied by anthers. 

In this paper we will propose a new finite element which is a bridge between con¬ 
forming and nonconforming hnite element.This new element piecewise quadratic 
and quasi conforming.We have studied error estimation for Dirichlet’s prob¬ 
lem,and observed that this element does not give optimal convergence rate,which 
has been generally considered as a major drawback against the use of this el¬ 
ement.Though the element is piecewise quadratic but it does not require two 
point continuity restriction on each interface of r^,which is needed for the above 
mentioned piecewise quadratic nonconforming element m .We shall show here 
that these elements are in fact very simple to use and they are nothing but 
combination of nonconforming element with incomplete P^ - conforming 
element. 


2 Discretization 

2.1 Continuous Problem 

We consider the following model problem 


—Au = / in D 
u = ud on 


( 1 ) 


where D C 
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Figure 1: NC1-C2 Element 
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2.2 Notation 

Let Th be a conforming triangulation of n.The subscript h refer to the maximum 
element size h = uvaxk^rh^k, where hk is the diameter of an element K £ Th-S^ 
is the set of the edges va Th, n is the unit outward normal along dK and the 
jump [u] across an edge e is a vector defined as follows- 

Let e be an interior edge shared by two triangles Ki and K 2 in , and Wj = 
u]\Kj for j=l,2 . We define on e 


[w] = wirii + W2n2 


where Uj is the unit normal of e pointing towards the outside of Kj .If e is an 
edge on the boundary of 12,then we define on e 


[w] = cun 


where n is the unit outer normal of e pointing towards the outside of 12. 

2.3 Weak Formulation 

Find u £V = Ftp (12) such that 


a{u,v) = F{v) y V GV = iFo(f^) 


where 



Q is convex polygon and / G L^(f2) therefore u G iF^(f2) by elliptic regularity 
theory [16]. 
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2.4 NC1-C2 methods 

In order to define a nonconforming space,we introduce some further notation. 

141 = {u € L^iQ) : v|/c is linear V K € Th, vis continuous at the mid points of the edges ofr/j} 
The above space is basically pi -nonconforming space. 

We define = span{(/)i, 02 , 03};where 0i = 0 ioF^^ and Fk is Affine mapping 
from AT to K [Fig 2]. 0i, 02 , 03 are basis fn on reference triangle K corresponding 
to the vertices 61 , 62 , ^3 respectively, which is defined by 

01 = (-1 -I- 2a; -I- 2j/)(-l + x + y) 

02 = (2x — l)x 
03 = {2y - l)y 

4>i> 1 < i < 3 is continuous along edge on each element.In this paper we use the 
following finite element space- 

Vh ■■= Vh ® Vh, Vh ■■= K e LH^) : e Dj,} 

Finite element space Vt consists of piecewise quadratic function which is 
discontinuous along edge of each triangle except at mid points of edges [Fig 1]. 

A typical element w € V/j is demonstrated below- 

w = -I- where G and 

since is discontinuous along edges except at the midpoints and is contin¬ 
uous along edge, w is discontinuous along edges except at the mid points. 

On interior edges jump of typical element is reduced to Pi polynomial .Let e be 
an interior edge which is shared by two triangle and and oj\K‘^ 

are restrictions of uj on and respectively. 

[w] = a;|Ar^ — u!\K‘^ 

= -P uj'^\K^) 

= 

= - uj^\K'^) 

This space contains the space of continuous piecewise-quadratic and space of 
nonconforming piecewise-linear,since 

Vh = V^+Pl 

where P^ is piecewise-quadratic conforming finite element space. 

2.5 Discretization 

The space Vh is not continuous and hence it is no longer in must 

modify the variational form a/i(.,.) in the discretized problem.We define the 
following bilinear form on Vh + F, .where F is subset of V and exact solution 
belongs to F- 
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2.6 Consistency 

Let u €V^t C V such that u satisfies weak formulation,i,e 

a{u,uj) = F{u;) V w G F (7) 

Then it is obvious that 

ahiu,u}) = Fh{uj) y u} €Vh ( 8 ) 

Hence the discrete bilinear form is consistent.lt also satisfies Galerkin’s orthog¬ 
onality condition since 


ah{u -Uh,u}) =0y u} eVh 


(9) 


2.7 Discrete stability 

The discrete bilinear form ah(-, ■) enjoys discrete stability on 14, if there is Csta > 
0 such that 

II 11 ^ O-hivhjtOh) 

Csta II Vfi II < II II 


ah 


{vh,uJh)= 'Y] / Vvh-VoJh 


KGTh 


considering Vh = tOh 


Vh 


= E / 

JK 


KGTh ' 


Hence 


Vh = 


= ahivh,Vh) 
ahivh,Vh) 

II Vh II 

ah{vh,ujh) 

< sup.,ec.\{0} II ^^11 


( 10 ) 

( 11 ) 


( 12 ) 


(13) 


Hence ah{-, ■) satisfies discrete stability condition.This implies the discrete bi¬ 
linear form 

ah{u,uj) = Fh{uj) y u,u} €Vh (14) 

is well posed i.e. the discrete bilinear form has unique solution. 


3 Apriori Error Estimation 

Lemma 1 Assume dim Vh < oo.Let ah{--) be a symmetric positive definite 
bilinear form on H -I- 14 which reduces to a(.,.) on V.Let u €V solve 

a{u, v) = F(v) V V € V 

where F € V' Cl Vj(.Let Uh G 14 solve 

ah(uh,v) = F{v) y V eV 
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Then 


II U ^ ' f II II I \^h{u Uli,U!) /I c:\ 

\\u-Uh |U< inf„ey^ II M - \\h +sup^gv^\{o} []—7]- (15) 

II ^ ll?i 

where || . ||?,= y^ahi-, ■) 

proof:- Let Uh € t4 satisfies 

ah{uh,v) = ah{u,v) V n € t4 


which implies that 


Then 


ah{uh - u,v) = 0\/ V eVh 

II M - Uh |U= inf„6Vfc II M - n ||?i 


\\U-Uh lU <11 M - Uh lU + II Uh - Uh II 
<11 u-Uh lU +sup^ev^\{o} 


\ah{uh - Uh,u})\ 

II ^lU 


Lemma 2 Let K be an arbitrary element of conforming triangulation r^i.Then 
the following inequality holds 

|er' II C Ili2(,)< II C fmK) +\C\miK)) V C e H\K) (16) 

where |e| denotes the length of edge e C dK, Hk =diam K,and the positive 
constant depends only on the chunkiness parameter of K. 


Proof: See the details in [T] 


Lemma 3 Let all assumptions of Lemma 2 hold and lo be an arbitrary 
element of 14 -Then 


kl II [w] ||i 2 (e)< C ^ (17) 

where [w] denotes jump of w along edge e € e^. 

Proof: Using lemma 2 we can write 

|e| ^ II M ||i2(g) - C" ^ (hj/ II Lo ||i2(;f) +l^lffi(_ftr)) 

KGTe 

|e| II H WUe) < C- E II ^ 11^4^) +\^\M^hHK)) 

KGTe 

Where is the set of triangles in Th containing e on their boundaries. Again 
[cli] =0 at midpoint of each edge of K hence first part of of right-hand side will 
be vanished.Therefore we have 

|e| II H lli^e) < C' E 

KGTe 

KGTc 


7 





Figure 2: Affine Mapping 
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Lemma 4 Let B be a ball in such that is star-shaped with respect 
to B and such that its radius p > {^)pmax- Let Q^u be the Taylor polynomial 
of order m of u averaged over B where u G W^(fl) and p > l.Then 

\u-Q"^u\wk(n) < Cm,n,'^d^~'^\u\w^in) k = 0,1, (18) 

where d=diam(fl) and Pmax=svcp { p : is star-shaped with respect to a ball 
of radius p } 

Proof: See the details in [T] 


3.1 

The important ingredient in the error analysis is a bound on the approximation 
error || u — uj || where uj € V/j is a suitable interpolation which agrees with u at 
mid point of each edges of of exact solution u.The interpolation operator is 
defined at the element level. We just require the local approximation property 

\u - ui\h-^{k) < Ch^j^^~"\u\HP+i(K) y K GTh ,s = 0,1,2 

It will be useful to define it explicitly.lt is defined in two steps. 

Let K be the reference triangle with vertices 6 i, 62 , ^3 whose coordinates are 
( 0 , 0 ),( 1 , 0 ),( 0 , 1 ) respectively and rhi be the midpoint of the side joining i and 
i-|-l(modulo 3) vertices. 

We define 

I^{u) = ■u(toi) 04 -I- u(m2)(()5 -I- u{lJl3)(j)Q 
P{u) = {u(bi) - /^(m)(Si))(^i -b {uih) - I^{u){b 2))^2 + (Hh) - 1 ^(m)(& 3))<^3 
Finally we define interpolation as 

i{u) = i^{u) + 

6 
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(19) 





where Li for i=l,2,3,4,5,6 continuous linear functional. 

Now we will show that P 2 {K) is unisolvent with respect to these functionals, 
i.e.for an arbitrary polynomial p G P 2 {K) Li(p) = 0 implies p = 0. 

Since p G P 2 {K) implies that p can be written as linear combination of basis of 
P 2 {k) .Then 

6 

P = '^C,(j>i. 

i=l 

C 4 = p{mi) = L4 {p) = 0 

Cs = p{m2) = kip) = 0 

Ce = pijhii) = Lq{p) = 0 
P = Ci^i + C2^2 + 

again 

Cl =p(Si) = Ti(:p) = 0 

C 2 = p{rh2) = L2 {p) = 0 
C 3 = p{rh3) = L3 {p) = 0 

=> p = 0 


Similarly it can be shown that for arbitrary p G P2{K), I{p) = p . Let 

{K, P2{K), S) be an affine finite element of {K, P2{k), S) where S = {Li, L2, L3, L4, L5, Lq} 

and S = {Li, 1/2, 1 / 3 ,1/4,1/5, Le} ■ 

Then for all v G we have 

II D^{u - hu) hHK)< C II \\l^(k) (20) 

where s < to + 1 and m=0,l,2. 


3.2 


We want to find out || u — Uh |U Using lemma 1 we can write 

II II ^ • r II II I (^h{u Uhjtu) 

\\u-Uh |U< mf„gy„ \\u-v \\h +sup^gy^\{0|--ji- 

II ^ lift 

\\u-v ||/i <11 u-ui ||/i 

< C k \u\h3{q) by ( 20 ) 

Now we have to estimate 


sup<^eVh\{o} 


|afe(u-«h,(/j)| 

II^^IU 


( 21 ) 
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ah{u — Uh,uj) = 2. / '^u.Vujdx— / fujdx 

KGTh 

= 22. [ / Vw-wncis — / AuwdO] — / fujdVl 
dk J K J K 


kGTh 


- E L Vu.wnds 


( 22 ) 


KGTh 


IdK 


= '^2 / VM.[a;]fis 


eg£'“ 


Again we have the following estimate 

E/ VM.[a;]iis = ^2 f - Ce)ne.[w]ds 


(23) 


eG£' 




e€£' 


h 


since /^[w] = 0 by (4). 

''22 / Vu.[w](is = '22 / (Vw.ne — Ce)ne.[w](is 

< XI 1 ®!”^ II Vu.ne-Ce ||L2(e) |e|i II M ||L2(e) 

< (22 |e|-' II Vu.rie-Ce ||i2(,))^(X |e| II H lli^(e))’ 


< ClX inmK(^rAh~2 II Vu.ne -Ce ||i 2 (T) +\ u \ h 2 i ^ k )) 2{22 X by (16,17; 

KGTe 


< Ch\u\H2(n) II uj ||/i 


Hence we have 


\ah{u - Uh,uj)\ <Ch |w|H 2 (n) || w ||/j 
\ah{u-Uh,uj)\ 




<Ch |u|// 2 (q) 


u-Uh iu< C h |m|h 2 (q) using (21) 


(24) 


3.3 L^-Error Estimate 

Let rj € iJ^(H) p| iJg (11) satisfy 


a{ri,v) = / v{u - Uh) V V € 
JVt 


and r]h G 14 satisfy 


ah{r]h,v)= / v(u-Uh)^v&Vh 
Jn 
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Again we have 

\\U-Uh 111,2(0) = J^{u-Uhf 

= {U- Uh){u - Uh) 

JQ 

= / u{u-Uh)- / Uh{u-Uh) 

Jo, Jci 

= a{u,r]) - ah{uh,'r]h) 

= ah{u,r]) - ah{uh,r]) + ah{uh,Tl) - ah{uh,r]h) 

= ah{u - Uh, v) + ah{uh, V - rih) 

= ah{u -Uh,r]- Vh) + ah{u - Uh,r]h) + ah{uh,ri - rjh) 
Using estimation (24) we can write 

ah{u-Uh,ri-r]h) ^Wu-UhWhWrj-VhWh (using Cauchy-Schwarz ) 

< C |u|/i-2(o) |?7|/i-2(Q) 

Where C is generic constant. 

ah{u - UhjVh) = ah{u - Uh,r]h - rji) + ah(u - Uh,r]h) 


( 25 ) 


( 26 ) 


ah 


ah{u - UhjTjh - rji) <11 u-Uh |U || ry/i “ Vi lU 

< C |m|h2(o) |r7|//2(o) 

II Vh - VI lU =\\Vh-V + V-Vi \\h 

<11 Vh-V \\h + II ?7-’7/ Ik 

< Cl h \v\H^(n) +C 2 h \v\H'^{n) 
= C h \v\H^(n) 

{u-Uh,rii) = '^ / VM.[r 7 /] ds by {22) 

eee'* 

= f Vu-lvi — v] ds 




= ^ / (Vw.ne - Ce) ne.[vi - v] ds 




Where Cg is arbitrary constant and we have used the following two identity 

W = o & /ek/] = o 

ah{u - Uh,Vl) < X! 1 ®!”^ II - Ce ||L2(g) |e |5 \\ v - Vi h^ie) 

eGe'* 

< C min/fgT-,.(/i^^infc,.eR || Vu.ng — Ce Ili2(;^;) +\u\'jh2(^K'f)]‘ 
^ lY. 1®!^ Y \\r]-Vi WhiK) +\V - Vi\m(K))]^ 

egjjit iCGTe 

< C |m|h 2 (o)(/i^ IvIh^ci)) using (18,20) 
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Therefore 


Again we have 


ah{u-Uh,r]h) < C h'^ \T]\H^(^n) 


ahiuh,r]- ilh) = ah{r] - r]h,Uh) 

= ahiv - r]h,Uh - ui) + ahiv - r]h,ui) 


(27) 


(28) 


ahiv - r]h,Uh -ui) < \\r]- Vh |U|| Uh - ui ||/i 


and 

II Uh - Ul \\h =\\ Uh - u + u - Ul \\h 

< \\ U - Uh \\h + \\ U - Ul ll/i 

< C h\u\H2(Q,) 

Therefore 

ahiv-Vh,Uh-ui) < C /i^|77|//2(o) |■u|//2(o) 

We can write second term of (28) as 

Uhiv -Vh,ui) = ^ / y-q-lui] 

= ^ / (Vr/.ne - Ce) ne.['u - Ul] ds 

where we have used the following two identity 

[■u] = 0 & [ui] = 0 ,e is an interior edge. 

Therefore 

Uhir] - ilh, ui) = '^ / (V? 7 .ne - Ce) ne.[u - uj] ds 

ede’’- 

^ X! lel"^ II V?7.ne - Ce ||L 2 (e) |e | ^ || u - || ^ 2 (e) 

< C |?7|//2(f2)(/l^ |u|ff2(Q)) 

Hence 

ahiuh,r]- Vh) = ahiv - r]h,Uh - uj) + ah{r] - r]h,ui) 

< C h'^ hlff2(n) |M|_f/2(f2) 

Using estimations (26),(27),(31) in (25) we can write 

\\u-Uh Ili 2 (n)< C Ivlmici) \u\m(n) 
using elliptic regularity we can write 

II ?? ||ff2(n)< C \\u — Uh\\L^(n) 

Therefore 

\\ u — Uh ||L 2 (n)< C |w|_f/ 2 (n) 


(29) 


(30) 


(31) 


(32) 

(33) 
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4 Numerical Test 


In this section we perform grid convergence studies for the proposed ’NCl- 
C2’ method.The new method was implemented using penalization technique.We 
estimate the experimental order of convergence by the formula 


Eoc = log{ 


E{2h) 

'W) 


) 


where E(h) =|| u — ut || is the error in the specified norm.The result indicate 
the same convergence behaviour for NC1-C2 method and and p^-nonconforming 
method. We have glued second order conforming element with first order noncon¬ 
forming element suitably and observed that first order nonconforming element 
dominated second order conforming element and we have obtained an EOC of 
2.0 in the norm and 1.0 in the norm which is same as nonconforming 
method.EOC stands for experimental order of convergence. 


Stationary Diflfusion Problem : We consider the Poisson equation 

- Am = / on O = (0,1) X (0,1) (34) 

with homogeneous boundary conditions and the right hand side / = 2t:‘^ sin('Kx)sin{'Ky). 
The exact solution is given by u(x,y) = sin{'Kx)sin(TTy). 
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